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ABSTRACT 

New kinematic data and modeling efforts in the past few years have substantially expanded and 
revised dynamical measurements of black hole masses (A/.) at the centers of nearby galaxies. Here 
we compile an updated sample of 72 black holes and their host galaxies, and present revised scaling 
relations between M. and stellar velocity dispersion (cr), F-band luminosity (L), and bulge stellar 
mass (Mbuigo), for different galaxy subsamples. Our best-fitting power law relations for the full galaxy 
sample are log 10 (M.) = 8.32 + 5.641og 10 (cr/200kms- 1 ), log 10 (M.) = 9.23 + 1.11 log 10 (L/10 n L ), 
and log 10 (M.) = 8.46 + 1.05 log 10 (Mb u i g e/10 11 M@). A log-quadratic fit to the M. — a relation with an 
additional term of fa [log 10 (fj/200kms _1 )] 2 gives fa = 1.68± 1.82 and does not decrease the intrinsic 
scatter in M # . Including 92 additional upper limits on M, does not change the slope of the M. — a 
relation. When the early- and late-type galaxies are fit separately, we obtain similar slopes of 5.20 and 
5.06 for the M, — a relation but significantly different intercepts - M m in early-type galaxies are about 
2 times higher than in late types at a given sigma. Within early-type galaxies, our fits to M,(a) give 
M, that is about 2 times higher in galaxies with central core profiles than those with central power- 
law profiles. Our M m — L and M. — Afbuigc relations for early-type galaxies are similar to those from 
earlier compilations, and core and power-law galaxies yield similar L- and Mb u ig -based predictions 
for M t . When the conventional quadrature method is used to determine the intrinsic scatter in M,, 
our dataset shows weak evidence for increased scatter at Mb u i ge < 10 11 Mq or Ly < 10 10 3 L@, while 
the scatter stays constant for 10 11 < Mbuigo < 10 12 3 M and 10 10 3 < L v < 10 n ' 5 L o . A Bayesian 
analysis indicates that a larger sample of M. measurements would be needed to detect any statistically 
significant trend in the scatter with galaxy properties. 

Subject headings: 



1. INTRODUCTION 

Empirical correlations between the masses, M,, of su- 
permassive black holes and different properties of their 
host galaxies have proliferated in the past decade. Power- 
law fits to these correlations provide efficient means to 
estimate M. in large samples of galaxies, or in individ- 
ual objects with insufficient data to measure M # from 
the dynamics of stars, gas, or masers. 

Correlations between black hole masses and numer- 
ous properties of their host galaxies have been explored 
in the literature. These include scaling rel ations be- 
tween M« and stellar velocity dispersion (e.g., |Ferrarese 
fc Merritt|| 2"000l ~ 



Gebhardt et al. 



|200i[ |Tremaine et al ||2002j |Wy the 2006a EfPu 



20001 JMerritt & Eer- 



rarese 

20081 iGultekin et al.||2009a[ hereafter G09; |Schulze fc 



Gebhardt||2011 McConnell et al. 2011a Graham et al 



201 1J [Beifiori et al. 



2012 hereafter B12) a nd between 



M. and" the stellar mass of the bulge (e.g., |Magorrian 



et~al~1[T998l |Marconi fe Hunt|[2003l |Haring fc Rix||2004| 



Hup009T[Sani et al.|[2011| B12j^ Various scaling rela- 
tions between M t and the photometric properties of the 
galax y have also been examined: bulge optical luminosity 
(e-g 



hardt 



Kormendy fe Richstone 1995| |Kormendy fe Geb- 



et al. 
Marc oni 

Graham|2007nmip009| |Sani et al.||20TTj j 



2001; G09; |Schulze fe Gebharo!t||Mll |Mc(%n5eII 
' 2011a[ B12V bulge near- infrared luminosity (e.g., 



fe Hunt||2003[ |McLure fc Dunlop|[2002 2004 
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ity (e .g., Kormendy & Gebhardt 2001 Kormendy et al. 



201 1| B12), and bulge concentration or Sersic index (eTg 



Graham et al. 2001 



a larger scale, corre 



Graham fc Driver 2007 B12). On 



ations between M, and the circular 



velocity or dynamical mass of the dark matter halo have 



been reported as well as disputed (e.g., Ferrarese 



Baeset al. 2003; Zasov et al. 2005 Kormendy & Bender 



2002 



2011[|Volonteri et al.|2011| B12)7More recently, M. has 
been found to correlate with the number and total mass 



of globular clusters in the host galaxy (e.g., Burkert fc 
Tremairie||20T0| |Harris fc Harri^|20TT| |Sadoun fc Colm 



2012[) . In ea rly-t ype galaxies with core profil es, Lauer 



et al.| ( |2007a[ ) and |Kormendy fc Bender] ( |2009[ ) have ex- 
plored correlations between M, and the core radius, or 
the total "light deficit" of the core relative to a Sersic 
profile. 

Recent kinematic data and modeling efforts have sub- 
stantially expanded the various samples used in all of 
the studies above. In this paper, we take advantage of 
these developments, presenting an updated compilation 
of 72 black holes and their host galaxies and providing 
new scaling relations. Our sample is a significant update 
from t wo recent compilations by G09 and |Graham et ah] 
( |2011| . Compared with the 49 objects in G09, 27 black 
holes in our present sample are new measurements, and 
18 masses are updated values from better data and/or 
more sop histicated mode l ing. Compared wit h the 64 
objec ts in Graham et al. (2011) (an update of Graham] 
20081, 35 of o ur black hole m asses are new or updates. 
The G09 and Graham (2008) samples differ by only a 
few galaxies, based on the authors' respective judgments 



about which dynamical measurements are reliable. The 
most significant upd ates in our sample are ga l axies with 
extremely high M. pien fc Gebhardt||2010l |Gebhardt| 

HKuslill 



terials to McConnell et al. (2011a 



et al. 2011 McConnell et al. 2011b 



and galaxies with some of the s mal 



2012 



poll 

est observed central 



black holes ([Greene et al.||2010| jNowak et al.||2010| |Kor^[ A), we have adopted the value M. 



mendy et al.] 2011 " 


Kuo et al. 






includes updated c 


Listances to 44 g 



)ur present sample 



The current sample 
includes on e new measurement of M, fro m |McConnell| 
et al. (2012), seven new measurements from Rush (2 012^7 
and t wo updated measurements (N GC 4594, |Jardel et al 
20TT| NGC 3998,|Walsh et al.|2012|). For NGC 5128 (Cen 

5.9±^ x 10 7 M w 



on surface brightness fluctuation meas urements ( Tonry 
eTaLpOOlj [Blakeslee et aLp009l [2010 I 



We focus on three frequently studied scaling relations: 
M. vs. stellar velocity dispersion (a), V^-band bulge lu- 
minosity (L), and stellar bulge mass (M bu i gc ). As re- 
ported below, our new compilation results in a signifi- 
cantly steeper power law for the M, — a relation than 
in G09 and the recent investigation by Beifiori et al. 
(2012; B12), who combined the previous sample of 49 
black holes f rom G09 with a large r sample of upper limits 
on M. from Beifiori et al.| ( 2009 ) . We still find a steeper 
power law than G09 or B12 when we include these up- 
per limits in our fit to the M, — a relation. We have 
performed a quadratic fit to M.(cr) and find a marginal 
amount of upward curvature, similar to previous investi- 
gations (Wyithe 2006a,b; G09). 

Another important measurable quantity is the intrin- 
sic or cosmic scatter in M, for fixed galaxy properties. 
Quantifying the scatter in M. is useful for identifying 
the tightest correlations from which to predict M. and 
for testing different scenarios of galaxy and black hole 
growth. In particular, models of stochastic black hole 
and galaxy growth via hierarchical merging predict de- 
creasing scatter in M« as galaxy mass increases (e.g., 
Peng||2007| |Jahnke fc Macci6||2011 ). Previous empirical 
studies of the black hole scaling relations have estimated 
the intrinsic scatter in M. as a single value for the entire 
sample. Herein, we take advantage of our larger sample 
to estimate the scatter as a function of <r, L, and Mb u i g c- 
In §2 we summarize our updated compilation of 72 
black hole mass measurements and 35 bulge masses from 
dynamical studies. In §3 we present fits to the M. — a, 
M. — L, and M. — Af^uigo relations and highlight subsam- 
ples that yield interesting variations in the best-fit power 
laws. In particular, we examine different cuts in a, L, 
and Mbuige, as well as cuts based on galaxies' morpholo- 
gies and surface brightness profiles. In §4 we discuss the 
scatter in M, and its dependence on er, L, and Mb„i gc . 
In §5 we discuss how our analysis of galaxy subsamplcs 
may be beneficial for various applications of the black 
hole scaling relations. 

Our full sample of black hole masses 
and galaxy properties is available online at 
http://blackhole.berkeley.edu. This database 
will be updated as new results are published. Investi- 
gators are encouraged to use this online database and 
inform us of updates. 



2. AN UPDATED BLACK HOLE AND GALAXY SAMPLE 

Our full sample of 72 black hole mass es and their host 
galaxy properties are listed in Table AT which appears at 
the end of this paper. The corresponding M, versus a, L, 
and Mbuigc are plotted in Figures TJ[3 This sample is an 
update of our previous compilation of 67 dynamical black 
hole measurements, presented in the supplementary ma- 



(at a distance of 4.1 Mpc) from Cappellari et al. (2009). 

We have removed three galaxies whose origina. 
surements ha ve exceptional complications 



mca- 



Lodato & 



Bertin (2003) measured non-Keplerian maser velocities 
in NGC 1068 and estimated M # by modeling a self- 
gravitating disk. Still, other physical processes might 
reprod uce the observed maser motions. [Atkinson et al. 
(2005) reported a measurement of M. in NGC 2748 but 
noted that heavy extinction could corrupt their attempt 
to locate th e center of the nuclear gas disk. IGebhardt 



et al. ( 2003 1 justified classifying the central point source 
of JN GC 7457 as an active galactic nucleus, but their ar- 
guments permit the central mass to be shared by an ac- 
creting black hole and a nuclear star cluster. 

Additionally, we have updated the distances to 44 
galaxies in our sample. For 41 galaxies, we adopt surface 



bright ness fluctuation measureme nts from Tonry et al. 
( |2001| and |Blakeslee et al.| j2009|), with the corrections 
suggested by jBlakeslee et ai.| l [Mo| ). For M31 and M32, 
we adopt the Cepheid variable distance of 0.73 Mpc from 

For NGC 4342 we a dopt the 
pc from Bogdan et al. (2012). Other 
measured quantities are scaled ac cord ingly: M. oc D, 

includes the up- 



Vilardell et al 
distance of 23 



M 



L cx D , ana ivibui, 
dated values for all quantities. The new galaxy distances 
and rescaled M. only have a small effect on our fits to the 
black hole scaling relations. For other galaxy distances 



and Mbuigc oc D. 



Table 



Al 



we assu me Hq = 70 km s Mpc , as in McConnell et al 



( |2011a| ). 

For the M, — a relation, we also consider upper lim- 
its for M, in 89 galaxies from B12, plus 3 new upper 



limits flSdiulze fc Gebhardt||2011[ |Gultekin et"aLl|2011 - 
McConn^ITTt™allI5nT5 1. Five additional galaxies in the 



B12 upper limit sample have recently obtained secure 
measurements of M, and are included in our 72-galaxy 
sample. As we discuss in §3, including upper limits re- 
sults in a lower normalization (intercept) for the M, — a 
relation but does not significantly alter the slope. 

For the M. — a relation, we consider two different defi- 
nitions of a. Both definitions use spatially resolved mea- 
surements of the line-of-sight velocity dispersion a{r) and 
radial velocity v(r), integrated out to one effective radius 
(r e ff): 



a 2 = 



_ /™ [a 2 (r) + « 2 (r)] I(r)dr 



frl 1 Wi- 



ll) 



where J(r) is the galaxy's one-dimensional stellar surface 
brightness profile. In G09 and most other studies, the 
lower integration limit r m i n is set to zero and sampled 
at the smallest scale allowed by the data. This defini- 
tion of cr, however, includes signal from within the black 
hole radius of influence, r in f = GM m a~ 2 . In some galax- 
ies, particularly the most massive ellipticals, a decreases 
substantially when spatially resolved data within r- m t are 
excluded. Setting r min = r in f produces an alternative 
definition of a that reflects the global structure of the 
galaxy and is less sensitive to angular resolution. We 
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Fig. 1. — The M% — a relation for our full sample of 72 galaxies listed in Table \KT\ and at http://blackhole.berkeley.edu Brightest 
cluster galaxies (BCGs) that are also the central galaxies of their clusters are plotted in green, other elliptical and SO galaxies arc plotted 
in red, and late-type spiral galaxies are plotted in blue. NGC 1316 is the most luminous galaxy in the Fornax cluster, but it lies at the 
cluster outskirts; the green symbol here labels the central galaxy NGC 1399. M87 lies near the center of the Virgo cluster, whereas NGC 
4472 (M49) lies ~ 1 Mpc to the south. The black- hole masses are measured using the dynamics of masers (triangles), stars (stars) or gas 
(circles). Error bars indicate 68% confidence intervals. For most of the maser galaxies, the error bars in Mm are smaller than the plotted 
symbol. The black dotted line shows the best-fitting power law for the entire sample: log 10 (M./ Mq) = 8.32 + 5.64 log 1 g(cr/200 kms - ). 
When early-type and late-type galaxies are fit separately, the resulting power laws are log lo (M./M ) = 8.39 + 5.20 log 10 (o-/200 kms" 1 ) 
for the early-type (red dashed line), and log 10 (Af./ Mq) = 8.07 + 5.06 log 10 (cr/200 kms -1 ) for the late-type (blue dot-dashed line). The 
plotted values of cr are derived using kinematic data over the radii r; n f < r < r c g. 



compare the two definitions of a for 12 galaxies whose 
kinematics within rj n f are notably different from kine- 
matics at larger radii. As shown in Table [ll excluding 
r < r in{ can reduce a by up to 10-15%. Ten of the 12 up- 
dated gala xies are mass ive (a > 250 km s -1 using either 
definition) . |Rusli| ( |2012[ ) presented seven new stellar dy- 
namical measurements of M. along with central velocity 
dispersions. We have used the long-slit kinematics from 
Rusli ( 2012[ ) and references therein to derive a according 



to Equation [T] our a values appear in Tables Al and [T] 
For the M, — Mbuige relation, we have compiled the 
bulge stellar masses for 35 early-type ga laxies. Among 
them, 13 bulge masses are taken from Haring & Rix 



(2004), who used spherical Jeans models to fit stellar 
kinematics. For 22 more galaxies, we multiply the V- 
band luminosity in Table |A1| with the bulge mass-to- 



light ratio (MjV) derived from kinematic s an d dynam- 
ical modeling of stars or gas (see Table Al for refer- 
ences). Where necessary, M/L is converted to y-band 
using galaxy colors. The values of M b u i ge are scaled to 
reflect the assumed distances in Table [All 

Most of the dynamical models behind our compiled val- 
ues of Mbuige have assumed that mass follows light. This 
assumption can be appropriate in the inner regions of 
galaxies, where dark matter does not contribute signifi- 
cantly to the total enclosed mass. Still, several measure- 
ments are based on kinematic data out to large radii. 
Furthermore, some galaxies exhibit contradictions be- 
tween the dynamical estimates of M/L and estima tes of 
M/L from stellar population synthesis models ( e.g., Cap- 
pellari et al. 120061 IConroy & van Dokkum|2012|). For this 



reason, we adopt a conservative approach and assign a 
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Fig. 2. — The M. — L relation for the 44 early-type galaxies with reliable measurements of the V-band bulge luminosity in our sample. The 
symbols are the same as in Figure 1. The black line represents the best-fitting power- law log lo (M./M ) = 9.23 + 1.11 log 10 (L„/10 n L s ). 



minimum error of 0.24 dex to each value of Afb u i ge . The vided by galaxy properties 
corresponding confidence interval (0.58 - 1.74) x Afbuige 
spans a factor of 3. 



To test how well our Mb u i ge values represent the stel- 
lar mass of each galaxy, we also have fit the M, — A/buigc 
relation using a sample of 18 galaxies for which Mbuige 
is computed from the stellar mass-to- light ratio, M+/L. 
Our stellar Mb u i gG sample comprises 13 galaxies for which 
M+/L is measured from dynamical models including dark 
matter, plus five galaxies for w hich M+/L is derived from 
stellar population models by |Conroy fc van Dokkurn] 
( |2012 ). This sample yields a slightly steeper slope of 
1.34±0.15 for the Af. — Mb u i ge relation, versus a slope of 
1.05 ±0.11 for our 35-galaxy dynamical Mbuigc sample. 
The stellar Mb u i g0 sample also has substantially lower 
scatter in M. (see Table [2]). 



3. BLACK HOLE SCALING RELATIONS AND FITS 

In this section we present results for the fits to black 
hole scaling relations for the full sample of dynamically 
measured M. listed in Tabic |A1| the full sample of M, 
plus 92 upper limits on AT,, and various subsamples di- 



3.1. Fitting methods 
Our power law fit to a given sample is defined in log 



space by an intercept a and slope 0: 



logio M » 



0log lo X. 



(2) 



where M, is in units of M e , and X = cr/200kms , 
L/10 n L©, or Afbuigc/lO^i© for the three scaling re- 
lations. We have also tested a log-quadratic fit for the 
M, — a relation: 



log 10 M. =a + f3 log 10 X + fa [log 10 X] 



(3) 



where X = er/200kms _1 . Res ults fo r the quadratic fit 
are discussed separately in Sec. 3.2.6 below. 

For the power-law scaling relations, we have com- 
pared three linear regression estimators: MPFITEXY, 
LINMIX_ERR, and BIVAR EM. M P FITEX Y is a least-squares 
estimator by |Williams et al.| (|20To| ). LINMIX_ERR is a 
Bayesian estimator by |Kelly| ( |2007l rBoth MPFITEXY and 
LINMIX_ERR consider measurement errors in two variables 
and include an intrinsic scatter term, eo, in log(M.). 
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Fig. 3. — The M. — Mt, u i ge relation for the 35 early-type galaxies with dynamical measurements of the bulge stellar mass in 
our sample. The symbols are the same as in Figure 1. The black line represents the best-fitting power-law log 10 (M./ Mq) = 
8.46 + 1.05 logioCMt^/lO 11 M @ ). 



LINMIX_ERR can be applied to galaxy samples with upper 
limits for M # . For the M. — a sample with upper limits, 
we also use the BIVAR EM algorithm in the ASURV soft- 
ware package by Lavalley et al. ( 1992), which implements 
the methods presented inllsobe et al.l (119861). The ASURV 



procedures do not consider measurement errors, and we 
use this method primarily for comparison with B12. All 
three algorithms are publicly available 3 . 

For each of the global scaling relations and galaxy 
subsamples, we obtain consistent fits from MPFITEXY 
and LINMIX.ERR, although LINMIX_ERR usually returns 
a slightly higher value of eo. Table [2] includes the global 
fitting results from both methods. In Table [2] we also in- 
clude results from LINMIX_ERR in cases where eo is poorly 
constrained by MPFITEXY. For the M. — a relation in- 
cluding upper limits, the BIVAR EM procedure returns a 

3 The IDL source code for MPFITEXY is available at 
http://purl.org/mike/mpfitexy . THE IDL source code 
for L INMIX_ERR and dependent scripts is available at 
http://idlastro.gsfc.nasa.gov/ftp/pro/math . ASURV is 
available at http://www2.astro.psu.edu/statcodes/asurv . 



lower intercept than LINMIX_ERR, but the slopes from 
the two methods are consistent within errors. Recently, 
Park et al. (2012) investigated the M, — a relation using 
tour linear regression estimators, including MPFITEXY and 
LINMIX_ERR. All four estimators yielded consistent fits to 
empirical data, and MPFITEXY and LINMIX_ERR behaved 
robustly for simulated data with large measurement er- 
rors in a. 

3.2. M. - a Relation 

Our fits to M.(cr) for the entire galaxy sample and 
various subsamples are plotted in Figures [T] and [4^i, and 
summarized in Table [2] 

3.2.1. Full Sample 

Our full sample of 72 galaxies yields an intercept a — 
8.32 ± 0.05 and slope /3 = 5.64 ± 0.32. When upper 
limits are added, the sample of 164 galaxies yields a = 
8.15 ±0.05 and = 5.58 ±0.30. The reduced intercept is 
a natural consequence of considering upper limits, while 
the slightly shallower slope is consistent within errors. 



5 



3.2.2. Early vs Late Type 

Fitting early- and late-type galaxies separately yields 
slightly shallower slopes: B = 5.20 ± 0.36 for early-types 
(red dashed line in Figure [T]) and ft = 5.06 ± 1.16 for late- 
types (blue dot-dashed line). The late-type galaxies have 
a significantly lower intercept: a — 8.39 ± 0.06 versus 
8.07 ± 0.21. Correspondingly, our fits predict M #iCar i y ~ 
2 M.jate at fixed a. Because most of the late-type bulges 
have low a, the split in intercepts leads to a steeper slope 
of 5.64 for the full sample. 



resolved) has the effect of decreasing a and increasing the 
slope of the M m — a relation. We have obtained new values 
of a for 12 galaxies in our sample by examining kinematic 
profiles from the literature and excluding data within 
r- m i; these galaxies are listed in Table [I] Ten of the 12 
galaxies have a > 250 km s -1 using either definition. 

In Table [2] we test how much the definition of a affects 
our fit to our full sample of 72 galaxies, as well as sub- 
samples dominated by massive early-type galaxies (see 



3.2.3. Core vs Power-law 

We also consider two subsamples of early-type galaxies 
classified by the slopes of their inner surface brigh tness 
profiles, 7 



dlogl/cilogr. Faber et al. (1997) and 



Lauer et al. (2007b) distinguished "power-law" galaxies 
with 7 > 0.5 from "core" galaxies with 7 < 0.3, although 
other studies have reported a continuous trend in 7 (e.g., 



Ferrarese et al. 2006 Glass et al. 2011). Core galaxies 



tend to be more massive and luminous than power-law 
galaxies, and there is some evidence that M, correlates 
with properties of the inner stel lar core (Lauer et al. 
2007a[|Kormendy fc Bender||2009[ ) . In our fits to M.{a), 
core galaxies have a significantly higher intercept than 
power-law galaxies (see Figure Ek): a = 8.53 ±0.11 ver- 
sus 8.24 ± 0.09. Our fits predict M. jCore ~ 2M,, p i at 
a ~ 200 kms , where the two populations overlap. The 
offset in intercepts plus the shallower slopes ((3 s» 4.5- 
4.8) for core and power-law galaxies combine to produce 
a steeper slope (/? ss 5.2) for the early-type M, — a rela- 
tion. 

3.2.4. Definitions of a 

As discussed in £[2j the value of a for each galaxy used 
in the M, — a relation depends on the spatial extent 
of the kinematic data. Excluding data within r in f (when 

TABLE 1 

Galaxies with multiple definitions of a 



Galaxy 


Ref. 




(0-r eS ) 


c (r inf -r e g) 






(") 


( kms -1 ) 


( kms -1 ) 


IC 1459 


1 


0.81 


340 


315 


NGC 1374 


2 


0.89 


203 


174 


NGC 1399 


3,4 


0.63 


337 


296 


NGC 1407 


5 


1.9 


283 


274 


NGC 1550 


6 


0.78 


302 


289 


NGC 3842 


7 


1.2 


275 


270 


NGC 3998 


8 


0.71 


286 


272 


NGC 4486 


9 


2.1 


375 


324 


NGC 4594 


10 


1.2 


240 


230 


NGC 4649 


11 


2.2 


385 


341 


NGC 4889 


7 


1.5 


360 


347 


NGC 7619 


12 


0.39 


324 


313 


NGC 7768 


7 


0.14 


265 


257 



Notes: References for kinematic data used to derive r- in { are (1) 
Cappellari et al. 2002; (2) D'Onofrio et al. 1995; (3) Graham et 
al. 1998; (4) Gebhardt et al. 2007; (5) Spolaor et al. 2008; (6) 
Simien & Prugniel 2000; (7) McConnell et al. 2012; (8) Walsh 
et al. 2012; (9) Gebhardt et al. 2011; (10) Jardel et al. 2011; 
(11) P inkney et al. 2003; (12) Pu et al. 2010. Although |Rusli| 
| |2012[ l report long-sl it kinematic measuremen ts for NGC 1374, the 
measurements from |D'Onofrio et ah] J 1995ft are m ore consistent 
with high-resolution SIN FU1N 1 data m |Kusli| ( |2012ft . 
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rows labeled "0-r e fr")- We find the slope of the global 
M, — a relation to change slightly from (3 — 5.64 ± 0.32 
in our fiducial sample (in which <7(r ln { — r e g) is used) 
to (3 = 5.48 ± 0.30 for the conventional definition of a 
with r nnn = in Equation ([I]). The latter is a fairer 
quantity to be compared with earlier studies, but the re- 
sulting M. — a relation is still significantly steeper than 
those reported in G09 and Bf 2. The definition of cr does 
not significantly affect our measurements of the intrinsic 
scatter in log(M.) (see Tableland §4). 

3.2.5. High vs Low a 

To search for possible systematic deviations of the 
Mm — a relation from a single power law, we divide the 
galaxies into low-cr and high-er subsamplcs, separated by 
a cutoff value cr C ut- We have tested numerous values of 
<7 cu t in search of robust trends. Our strongest finding 
is that the relation for the highcr-cr sample appears to 
steepen drastically when <7 cut 270 kms -1 . As shown 
in Table [2] (for cr cut = 275 kms -1 ), the MPFITEXY and 
LINMIX_ERR procedures both return nearly vertical rela- 
tions, with very large uncertainties in (3. This suggests 
a breakdown of the M. — a correlation as the galaxy 
population "saturates" at cr ~ 350 kms -1 . Saturation 
of the L — a and M. — cr relations has been predicted 
from observa tions and simulations of the most massive 
galaxies (e.g.,[B oylan- Kalchin et al.|2006| |Bernardi et al 



2007 



Lauer et al.]J2007a) 



The MPFITEXY procedure re- 
turns zero intrinsic scatter when fitting galaxies with 
a > 275 kms -1 ; this is an artifact of the large slope. 

One might suspect that the saturation in M, (a) results 
from the lower cr values we obtained for 12 galaxies after 
excluding data within r; n f. We have also fit the high-er 
and low-cr subsamples using the conventional definition 
of cr. We still find that the highest-er galaxies follow a 
very steep M. — cr relation, although this trend begins to 
appear at slightly higher values of cr cu t (exemplified by 
cr cut = 290 kms -1 in TableJJ. 

A weaker trend occurs for cr cut in the range 175-225 
kms . Here, the higher-er sample exhibits a steeper 
slope [fi ~ 6-7) than the lower-cr sample ((3 ~ 5). Still, 
the differences between the fits for the two subsamplcs 
are within ler error bars as the uncertainties in a and j3 
are large. This trend vanishes when we adopt the conven- 
tional definition of a (data from r = 0- r e s). As an exam- 
ple, Table|2]lists the fitting results for cr cut = 200 kms -1 . 

It is tempting to fit M.(cr) for narrow intervals in a 
(e.g., 150 kms -1 < a < 200 kms -1 ), but the intrinsic 
scatter in the M. — a relation drives these samples toward 
an uncorrelated distribution. For the current sample of 
Mm, a long baseline in a ( 100 kms -1 ) is therefore 
needed to determine the slope of the M, — cr relation. 

We have also tried fitting M.(cr) for subsamples de- 
fined by cuts in L and A/b u i ge ; two examples are listed in 
Table [2j Some of the high-L subsamples exhibit a steep 
Mm — c slope, but again with large uncertainties. 

3.2.6. A Log- quadratic Fit to M % — cr 

In light of evidence that the M. — a relation steepens 
toward high galaxy masses, we have also attempted to fit 
Mm (cr) as a log-quadratic function in Equation The 
coefficients a, /3, f3i, and intrinsic scatter e for our 72- 
galaxy sample are determined from a brute-force least- 
squares estimator similar to MPFITEXY. We find that the 



best-fit parameters are a = 8.28 ± 0.07, (3 = 5.76 ± 0.34, 
/3 2 = 1.68 ± 1.82, and eo = 0.38. Uncertainties in a, (3, 
and P2 are determined by assessing the one-dimensional 
likelihood function after marginalizing y 2 with respect 
to the other two parameters. We find P2 > with 82% 
confidence, slightly below the lcr threshold for a one- 
sided confidence interval. 

As suggested by the highly uncertain power-law slopes 
for high-cr galaxies, our measurement of upward curva- 
ture in Mm(o~) is marginal. Adopting a quadratic relation 
does not decrease the intrinsic scatter in M,. Our full 
sample updates the investigations of Wyithe (2006a, b; 
31 galaxies) and G09 (49 galaxies), who reported simi- 
lar confidence levels for a non-zero quadratic term. At 
the extreme end of the local galaxy velocity dispersion 
function (cr ~ 400 kms -1 ), our best quadratic fit predicts 
black hole masses ~ 40% higher than the best power-law 
fit. 

3.3. Mm — L and M, — Mbuigc Relations 

In Table |A1| we present V— band luminosities for 44 
galaxies and dynamically measured bulge masses for 35 
galaxies. For several of the late-type galaxies in our sam- 
ple, the literature contains one or more estimates of the 
bulge-to-total light ratio. Rather than judging between 
the various estimates, we present results for the early- 
types only. Our best fit values are (3 — 1.11 ± 0.13 
and a = 9.23 ± 0.10 for the M, — L relation, and 
/3 = 1.05 ± 0.11 and a = 8.46 ± 0.08 for the M. - M bu i ge 
relation. 

Additional fits to subsamples of these galaxies are 
listed in Table [2j The M, - L and M, - M bu i ge relations 
do not show statistically significant differences between 
core and power-law galaxies (see also Figures EJd and [4p). 

Figures [2] and [3] show that our M, — L and M, — Mb u i gc 
samples both appear to have a central knot, where black 
holes with 10 8 M G < M. < 10 9 M exhibit relatively 
weak correlation with L or -Mbuigc- This feature makes it 
difficult to interpret the fits to high-L and low-L (or high- 
-Mbuigo and low-Mb u i g e) subsamples. We find tentative 
evidence that the most luminous and massive galaxies 
(L > 10 10 - 8 L ; Mbuigc > 10 115 M Q ) have steeper slopes 
in Mm{L) and M # (Mb u i g e), as exemplified in Table 
Both samples are sparsely populated at the low-M. eh 

3.4. Comparison to Previous Studies 

The slope of the M, — cr relation reported in prior stud- 
ies has wavered between ~ 4 (e.g., Gebhardt et al. 2000 
Tremaine et al. 2002; G09; B12) and ~ 5 (e.g.,|Ferrarese 
k Merritt|2000l|Merritt k. Ferrarese|2001||Graham et al. 



2011 1. Our best-fit slope for the global M. — cr relation 
falls at the steep end of this distribution, while various 
subsamples exhibit a wider range of slopes (/3 ~ 3.8 to 
/3 > 12). In particular, the M # — a relation for our full 
sample is significantly steeper than those reported in G09 
(13 = 4.24 ±0.41) and B12 (j3 = 4.42 ±0.30). This steep- 
ening has occurred because the newest measurements of 
Mm in early-type galaxies (higher cr) mostly fall above the 
global Mm — cr relation, and the newest measurements of 
Mm in late-type galaxies (lower cr) mostly fall below the 
global relation. In addition to the significant discrepancy 
between the two subsamples' best-fit intercepts, both the 
early- and late-type M, — cr relations have steepened. 
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TABLE 2 

Power-law fits to black hole correlations 



Sample 




Method 


a 




£0 


M, — a relation 












All galaxies 


72 


MPFITEXY 


8.32 ± 0.05 


5.64 ±0.32 


0.38 


All galaxies 


72 


LINMIX.ERR 


8.31 ± 0.06 


5.67 ± 0.33 


0.40 ± 0.04 


All + upper limits 


164 


A C T TT» T r 


"7 70 _L n 1 o 


O.O / ± U.oz 




All + upper limits 


104 


T T HTM TV UDD 

LlWMlA_hriK 


o ik i n nc; 
O.10 ± U.UO 


O.Oo ± U.oU 


n a q J_ n f\A 
U.4o ± U.U4 


All galaxies (0-r e ff) 


72 


MPFITEXY 


8.29 ±0.05 


5.48 ± 0.30 


0.37 


G09 data (0-r off ) 


49 


MPFITEXY 


8.19 ±0.06 


4.12 ±0.38 


0.39 


Early-type 


Oo 


MPr 1 IhXY 


Q Qfl J_ n Oft 

o.oy ± U.Uo 


on J_ n Qft 
O.zU ± U.oo 


n q /i 
U.o4 


Early-type (0-r off ) 


53 


MPFITEXY 


8.36 ±0.05 


5.05 ±0.34 


0.33 


Late- type 


19 


MPFITEXY 


8.07 ±0.21 


5.06 ± 1.16 


0.46 


Power-law 


18 


MPf 1 IhXY 


q O/i _L n nn 

o.z4 ± u.uy 


A C 1 _l_ n 7Q 

4.0l ± U. to 


0.34 


Core 


OQ 

zo 


MPr 1 IhXY 


Q ctq 1 n 1 1 
o.Oo ± U.ll 


A 70 _L n 1 A 

4. r y ± U. ( 4 


U.oO 


Core (0-r eff ) 


28 


MPFITEXY 


8.50 ±0.11 


4.63 ± 0.68 


0.34 


u <~ zUU km s 


35 


MPf 1 IhXY 


o.oo ± U.lo 


f: i not; 
O.OO ± U.oo 


0.43 


\ onn i,-, ~ — 1 

(7 > zUU km s 


37 


MPf 1 IhXY 


o i a _i_ n 1 o 
o.lo ± U.lo 


a Ta _L n n 1 

O. it ± u.yi 


0.34 


a > 200 km s" 1 (0-r cff ) 


38 


MPFITEXY 


8.26 ±0.12 


5.70 ±0.78 


0.35 


<r <^ z / o km s 


00 


T T1IMTV T7DD 


O.OO ± U.U f 


0. it ± U.Ol 


U.4o ± U.Uo 


(T > z ( o km s 


1 7 


T TTlTlkJTTV TTT)r> 

LlNMlA_hRK 


7 nn _L o a o 
f .UU ± 2.42 


1 O Q _l_ 1 O C 

12. o ± 12. D 


q A _l_ n 1 1 
U.34 ± U.ll 


cr > 275kms~ 1 


17 


MPFITEXY 


2.47 ±3.17 


35.8 ± 16.5 


N/A 


C7 > 290 km s" 1 (0-r cff ) 


15 


LINMIXJLRR 


7.68 ± 1.26 


7.93 ±5.79 


0.32 ±0.11 


L < 10 lue Lq 


25 


MPFITEXY 


8.37 ± 0.08 


4.76 ± 0.55 


0.33 


t -i nl fl 8 r 

L > 10 1U ' 8 Lq 


19 


MPFITEXY 


8.13 ± 0.40 


7.19 ± 2.25 


0.37 


L > 1O 1O - 8 L (0-r off ) 


19 


MPFITEXY 


8.29 ±0.33 


5.83 ± 1.75 


0.36 


M bu i gc < 10 11 ' 5 M Q 


21 


MPFITEXY 


8.40 ± 0.09 


5.08 ± 0.70 


0.34 


M bu l ge > 10 11 - 5 Mq 


14 


MPFITEXY 


8.52 ± 0.47 


4.69 ±2.69 


0.46 


M bulgc > 10 n - 5 M Q (0-r cff ) 


14 


MPFITEXY 


8.61 ± 0.40 


3.80 ±2.09 


0.45 


M, — L relation 












Early-type galaxies 


44 


MPFITEXY 


9.23 ±0.10 


1.11 ± 0.13 


0.49 


Early-type galaxies 


A A 
44 


T T1IMTV rnn 


fl OQ 1 r> 1 r\ 

y.zo ± U.1U 


1 1 1 i n 1 \ 
1.11 ± U.14 


n c;o J_ n nft 
U.Oz ± U.Uo 


G09 data (early-type) 


32 


MPFITEXY 


9.01 ± 0.10 


1.17 ± 0.12 


0.36 


Power-law 


12 


MPFITEXY 


9.36 ± 0.72 


1.19 ±0.67 


0.68 


Core 


27 


MPFITEXY 


9.28 ± 0.09 


1.17 ± 0.22 


0.39 


L< 10 10 - 8 L o 


25 


MPFITEXY 


9.10 ±0.23 


0.98 ± 0.20 


0.54 


L > 1O 1O 8 L 


19 


MPFITEXY 


9.27 ± 0.13 


1.12 ±0.82 


0.47 


M bu i gc < 10 115 M Q 


18 


MPFITEXY 


9.25 ± 0.24 


1.13 ±0.23 


0.47 


M bulgc > 10 115 Mq 


13 


MPFITEXY 


9.24 ±0.10 


2.49 ±0.67 


0.30 


M, — M bulge relation 












Dynamical masses 


35 


MPFITEXY 


8.46 ± 0.08 


1.05 ±0.11 


0.34 


Dynamical masses 


35 


LINMIX.ERR 


8.46 ± 0.09 


1.05 ±0.12 


0.36 ± 0.08 


Stellar masses 


lo 


MDfTTUYV 

MPr 1 IhXY 


O.OO ± U.1U 


i Q/i i n 1 ^ 
l.o4 ± U.10 


U.l / 


JT UWtl law 


12 


MPFTTFXY 
iir r ± i r_jA i 




n 04 _|_ n qq 


0.50 


Core 


20 


MPFITEXY 


8.45 ± 0.15 


1.09 ±0.20 


0.28 


L< 1O 1O8 L 


19 


LINMIXJLRR 


8.44 ±0.14 


1.05 ±0.26 


0.45 ±0.13 


L> 1O 1O8 L 


12 


LINMIX.ERR 


7.66 ± 1.60 


1.92 ± 1.72 


0.38 ±0.19 


L> 10 1() - 8 L Q 


12 


MPFITEXY 


6.92 ± 1.05 


2.72 ± 1.12 


N/A 


M bulgc < 10 115 Mq 


21 


LINMIXJLRR 


8.54 ±0.15 


1.11 ± 0.28 


0.47 ± 0.12 


M bulgc > 10 115 Mq 


14 


LINMIXJLRR 


7.28 ± 1.19 


2.26 ± 1.33 


0.30 ±0.17 


M bulgc > 10 115 M Q 


14 


MPFITEXY 


7.03 ± 0.78 


2.53 ±0.85 


N/A 



Notes: For the M. — c relation, we fit log 10 (M.) = a + /3 log 10 (<r/200 km s _1 ). Subsamples designated (0-r e ff) define a using kinematic 
data over the interval < r < r c g . For all other subsamples, we define a using data over the interval rj n f < r < r c g. For the M. — L 
relation, we fit log 1 g(M.) = a + /3 log 10 (L/10 11 Lq). Luminosities are in V-band. For the M. — M bu i ge relation, we fit log 1 Q(M.) = 
a + /3 log 1 Q(M bu i ge /10 11 Mq). All fits except for the "stellar masses" line use the sample of bulges with dynamical masses. 
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Our fit to early-type galaxies is significantly steeper 
than the early-type fit by G09 ((3 = 5.20 ± 0.36, versus 
/3 = 3.96 ±0.42). This difference largely is due to several 
updates to the high-er galaxy sample: new measurements 
of M. ~ 10 10 M Q in the brightest cluster gal axies NGC 
4889 and NGC 3842 (jMcConnell et al.|2011a[ ), new mea- 



surements of M« > 10" Mf) in seven more galaxies (Rush 



et al. 2011 Rush 2012), and updated measurements in- 



creasing~M # in M87 and M60 ( |Gebhardt fc Thomas 2009 
Shen fc Gebhardt|2010[ ). Defining ct to exclude the black 
hole radius of influence further steepens the early-type 
galaxy sample by a small amount. I f we exclude th e re- 
cent additions by IMcConnell et al.l ( |2011a|b [ |2012[ ) and 
Rusli| ( |2012[ ), we obtain (3 = 4.77 ±0.36 for 42 early-type 
galaxies. Removing M87 and M60 further reduces (3 to 
4.55±0.37; in addition to revised black hole masses, these 
two galaxies exhibit some of the largest differences in ct 
in Tabled] 

Our fit for late-type galaxies is slightly steeper than 
G09 (/3 = 5.06 ± 1.16, versus j3 = 4.58 ± 1.58). This 
arises primarily from our exclusion of NGC 1068 and 
NGC 2748. 

O ur earlier compilation of a similar sample of 67 galax- 
ies ( jMcConnell et al.||2011a| ) gave a = 8.28 ± 0.06 and 
f3 = 5.13 ± 0.34 for the M. - a relation. Our M.(a) fit 
to the present sample of 72 galaxies has a steeper slope 
of 5.64 ± 0.32, largely due to the exclusion of NGC 7457, 
which had the lowest velocity dispersion (ct — 67kms~ 1 ) 
of all galaxies in the previous sample (see § [2] and Geb- 
hardt et al. 2003 for discussion of this galaxy's central 
massive object). If we include NGC 7457 in our present 
sample, we obtain a = 8.33 ± 0.05, (3 = 5.42 ± 0.31, and 
eo = 0.40, closer to our earlier results. 

Our M # — L and M, — Mb u i ge slopes are consistent with 
a number of previous investigations, including multiple 
bandpasses for L (e.g., Marconi & Hunt 2003; Haring 
& Rix 2004; McLure & Dunlop 2004; G09; Schulzc & 
Gebhardt 2011). For the M. - L relation 
(2011 ) report different M. — L 3 , 6lim and M. 



Sani et al. 
- Ly slopes 

and suggest that color corrections and extinction may be 
responsible for the difference. Their M # — £3.6^111 slope is 
0.93±0.10, while their M, — Ly slope ranges from 1.11 to 
1.40 depending on the regression method. These slopes 
are consistent with our M. — Ly slope of 1.11 ± 0.13. 

For the M, — Mb u i go relation, the latest compilation 
of 46 galaxies by B12 gives a slope of 0.79 ± 0.26. We 
note that their Mb u i ge values are virial estimates based 
on the galaxies' ct and r ff via Mb u i gc = 5. Oct 2 r e g/G. 
In comparison, our Mb u i gc values use the mass-to-light 
ratios obtained from dynam ical models. 

( |2012[ ) examined the M. — ct and 



Recently, Graham 



Mm — Mbui ge relations with separate fits to core and non- 



core galax ies, based on the galaxy sample of Haring <k 



Rix (2004) and updated black hole masses from Graham 
et al.f j 2011 ) . The non-core galaxies were found to follow 
a very steep M, — Mb u i ge relation (/? ~ 2), and there 
was virtually no difference in the M. — ct relations for 
core versus non-core galaxies. Our relative trends for 
core a nd power-law galaxies differ from those in | Graham] 
(2012). This is likely due to differences in the galaxy 
samples: our core galaxies include eleven galaxies with 
Mm > 10 9 Mp, that are absent from the sample used by 
Graham (2012). Our photometric classification of galax- 



ies also differs from Graham (2012). In particular, we 
classify the high-M bu i ge object NGC 6251 as a power- 
law galaxy, based on t he surface brightness profile of 
Ferrarese & Ford (1999). Excluding NGC 6251, we mea- 
sure f3 « 1.6 for power-law galaxies on the M. — -Mbuigc 
relation. 

4. SCATTER IN BLACK HOLE MASS 

For a given black hole scaling relation, the differences 
between the measured values of M m and the mean power- 
law relation are conventionally interpreted as a combi- 
nation of measurement errors and intrinsic scatter. We 
assume the scatter in M, to be log normal, and define 
the intrinsic scatter term eo such that 



X 



E 



[log 10 {Mm,i) -a- fai 



-M,i 



+ /3 2 



(4) 



where x = log 10 (CT/200kms _1 ) for the M, — ct relation, 
x = log 10 (L/10 Lq) for the M. — L relation, and 
x = log 10 (Mbui gc /10 Mq) for the M,-Mb u i gc relation. 
Here, 6m is the 1-ct error in log 10 (M.), and e x is the 1- 
ct error in x. For a given sample and power-law fit, we 
adopt the value of e for which xl — 1 (x 2 = ^dof)- G09 
tested several forms of intrinsic scatter in M, and found 
log normal scatter to be an appropriate description. 

Fitting the full galaxy sample for each scaling relation, 
we find the intrinsic scatter in log 10 (Af») to be eo = 0.38 
for the M.-ct relation, eo = 0.49 for the M t -L relation, 
and eo = 0.34 for the M.-Mb u i go relation (or eo = 0.17 
for Mm versus stellar Mb u i gc ). While it is tempting to 
conclude that Mb u i g o is the superior predictor of M. , the 
relative errors in Mb u i gc , cr, and L demand a more cau- 
tious interpretation. As noted in ^2j we have assumed 
that all Mbui ge values have an error of at least 0.24 dex. 
We have repeated our fits to M. (ATbui gc ) with a minimum 
error of only 0.09 dex. Fitting the full Afbuigc sample 
with this reduced error in Mb u i ge , we obtain a larger in- 
trinsic scatter (e = 0.39) as expected from Equation |4j 
while the slope and intercept of the fit do not change 
significantly. Similarly, our measurements of eo for the 
Mm — ct and M. — L relations depend in part upon the as- 
sumed errors in ct (> 5%, following G09) an d L (typically 
< 0.05 dex). In addition to evaluating eo, 
(20061 used earlier datasets to assess whia 



Novak et al 
ti correlation 

yielded the lowest predictive uncertainty in M., given a 
set of host galaxy properties with measurement errors. 
They found the AI t — ct relation to be marginally favor- 
able for predicting black holes with M, ~ 1O 8 M , but 
noted that uncertainties in the relations' slopes compli- 
cated predictions near the extrema of the relations. Our 
global Mm — ct relation has a steeper slope ((3 = 5.64) 



than the samples evaluated by Novak et al. ( 2006 1 , with 
/3 from 3.69 to 4.59. 

Beyond the intrinsic scatter e in black hole mass for 
the full sample, the dependence of eo on ct, L, and Mb u i gc 
is a useful quantity for constraining theoretical models 
of black hole assembly. Successive mergers are predicted 
to drive galaxies toward the mean M, — Mb u i go rela- 
tion, especially when black hole gro wth is dominated by 
black hole-black hole mergers (e.g., Peng 2007 Jahnke 
& Maccio 20111). Understanding intrinsic scatter in M, 
is also crucial tor estimating the mass function of black 
holes, starting from the luminosity function or velocity 
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Fig. 5. — Scatter in logjg M, for different intervals in cr, L, and -Mb u i gc . Black squares represent the intrinsic scatter, eo required 
to obtain x 2 = -^gal between the subsample of galaxies and the global scaling relation. Red diamonds represent the root-mean-squarcd 
residual for each interval, between log(M. ) and the global scaling relation. Blue stars with vertical error bars represent the Bayesian 
estimates for eo, obtained while fitting a separate scaling relation for each interval, (a) Scatter with respect to the Mm — cr relation, 
log 10 M m = 8.32 + 5.64 log 1( )(o-/200kms — 1 ), defining cr with data from r; n f to r cn -. (b) Scatter with respect to the M. — cr relation, 
log 10 Mm = 8.29 + 5.48 log 1 Q(cr/200 kms -1 ), defining cr with data from to r c ff. (c) Scatter with respect to the M. — L relation, log 10 M. = 
9.23 + 1.11 log^L/lO 11 Lq). (d) Scatter with respect to the Mm - Mbulge relation, log 10 Mm = 8.46 + 1.05 log 10 (A^ bulgc /10 11 M ). 



dispersion function of galaxies (e.g., Lauer et al. 2007c; 
G09). 

Figure [5] illustrates how eo varies across each of the 
Mm — cr, Mm — L, and M. — Mb u i g c relations for our up- 
dated sample of measurements. For each relation, we 
construct three or four bins containing equal numbers 
of galaxies and perform multiple estimates of the scat- 
ter in each bin. One estimate, represented by blue stars 
in Figure [5] is to perform an independent Bayesian fit 
(with LINMIX_ERR) for the scaling relation in each bin, 
and assess the posterior distribution of eo- This method 
provides uncertainties for eo in each bin, but the fits to 
narrow data intervals typically yield poor estimates for 
the slope and intercept. A second estimate, represented 
by black squares in Figure [5] is to compute eo as in Equa- 
tion^ using the global fit to the scaling relation to define 
the same values of a and /3 are for all bins. The third and 
simplest estimate, represented by red diamonds in Fig- 
ure^] is to evaluate the root-mean-squared (RMS) resid- 
ual between log(Af.) and the global scaling relation. The 
eo term in Equation [4] provides a reliable assessment of 
intrinsic scatter only if random measurement errors are 
small: measurements with large uncertainties can yield 
xt < 1 with no intrinsic scatter term. In comparison, 
the RMS estimate has no explicit dependence on mea- 
surement errors. 



Figures [5^ and [5}d illustrate the scatter in M. as a 
function of a. Two definitions of a for the 12 galaxies 
in Table [l] are shown for comparison, where a is com- 
puted from data between radii r- ln f and r e s, or between 
and r c ff . Considering the large error bars in eo from the 
Bayesian fits, we find no significant variation in M, with 
respect to a. The high- and low-mass ends of the M. — a 
relation both exhibit ~ 0.3-0.4 dex of scatter, regardless 
of how a is defined or how scatter is estimated. 

For the M. — L and M, — Mbulge relations, we find 
possible evidence that galaxies with low spheroid lu- 
minosities (L < 10 103 L Q ) and small stellar masses 
(Mbulge < 10 11 M Q ) exhibit increased intrinsic scatter 
in Mm- Yet the scatter appears constant for galaxies 
above this range, spanning 10 10 3 L Q < L < 10 115 L Q 
and lO n M < M bu i g c < 10 12 - 3 M Q . More measure- 
ments in the range Mb u i gc ~ 10 8 — 10 10 M are needed 
to reveal whether intrinsic scatter in M, varies system- 
atically across an extended range of bulge luminosities or 
masses. The Bayesian estimates of eo in each bin have 
large uncertainties; adopting this method, we do not de- 
tect a significant change in scatter for any interval in cr, 
L, or Mbuigc- 

The Local Group galaxy M32 is separated from the 
other early-type galaxies in our sample by almost an or- 
der of magnitude in L, and more than an order of mag- 
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nitude in Mbuige- We have excluded its contribution in 
Figures [HJ: and [5ji, so that the sizes of the leftmost bins 
better reflect the sampled distributions of L and Mbuige- 
Including M32 does not substantially change the amount 
of scatter in the lowest- L and lowest-Mb u i ge bins. 

Although our three estimates of scatter do not yield 
the same absolute values, their qualitative trends as a 
function of a, L, or -Mbuige are very similar. The simi- 
lar behavior of eo and RMS indicates that variations in 
measurement errors are not responsible for the apparent 
trends in intrinsic scatter. 

5. SUMMARY AND DISCUSSION 
We have compiled an updated sample of 72 black 



hole masses and host galaxy prop ertie s in Table AT 
a more detailed version of Table I All is available at 



http://blackhole.berkeley.edu. Compared with the 
49 objects in G09, 27 black holes in our sample of 72 are 
new measurements and 18 masses are updates of pre- 
vious values from improved data and/or modeling. Our 
present sample includes updated distances to 44 galaxies. 

We have presented revised fits for the M. — cr, M. — L, 
and M. — Mb u i ge relations of our updated sample (Ta- 
ble]^ and Figures T|3). Each relation is fit as a power- 
law:log 10 (Af.) = a + /3 log 10 (A). Our best fit to the full 
sample of 72 galaxies with velocity dispersion measure- 
ments (X = cr/200kms" 1 ) is a = 8.32 ± 0.05 and /3 = 
5.64 ± 0.31. A quadratic fit to the M, — a relation with 
an additional term /3 2 [log 10 (A)] 2 gives (3 2 = 1-68 ± 1.82 
and does not decrease the intrinsic scatter in M # . In- 
cluding 92 additional upper limits for M, decreases the 
intercept but does not change the slope: a = 8.15 ± 0.05 
and (3 = 5.58 ± 0.30. 

For the 44 early-type galaxies with reliable T^-band 
luminosity measurements (X = Ly /10 11 £©), we find 
a = 9.23±0.10 and /3 = 1.11±0.13. For the 35 early-type 
galaxies with dynamical measurements of the bulge stel- 
lar mass (X = Mbuigc/10 11 M ), we find a = 8.46 ±0.08 
and P = 1.05 ±0.11. 

We have also examined the black hole scaling rela- 
tions for different subsamples of galaxies. When the 
galaxies are separated into early- and late-types and fit 
individually for the M t — a relation, we find similar 
slopes of (3 = 5.20 ± 0.36 (early-types) and 5.06 ± 1.16 
(late- types) . The intercepts, however, differ significantly: 
a = 8.39 ±0.06 for the early types, a factor of ~ 2 higher 
than a = 8.07 ± 0.21 for the late types. The steep global 
slope of 5.64 is therefore largely an effect of combining 
different galaxy types, each of which obey a shallower 
M. — cr relation and different intercepts. 

When the early-type galaxies are further divided into 
two subsamples based on their inner surface brightness 
profiles, the resulting Af, — a relation has a significantly 
larger intercept for the core galaxies than the power- 
law galaxies (Table [2] and Figure |4^l) . The slopes of the 
M. — L and M. — Aibuigc relations do not show statisti- 
cally significant differences between core and power-law 
galaxies, but M. follows L and -Mbuige more tightly in 
core galaxies than power-law galaxies (Table [2]). 

In the literature, the exact value of the M7~— a slope 
has been much debated by observers and regarded by 
theorists as a key discriminator for models of the assem- 
bly and growth of supermassive black holes and their host 
galaxies. We suggest that the individual observed M, — a 



relations for the early- and late- type galaxies provide 
more meaningful constraints on theoretical models than 
the global relation. After all, these two types of galax- 
ies are formed via different processes. Superficially, our 
measurement of /? = 5. 64 ±0.32 for the global M. — a re- 
lation favors thermally-driven w ind models that predict 
/3 ~ 5 (e.g., Silk & Rees 1998 ) over momen tum-driven 



wind models with [3 ~ 4 (e.g., Fabian 19991. However 



the intercept of the empirical Ai # — a relation is sub- 
stantially higher than int ercepts derived from thermally- 
driven wind models (e.g., |King||2010a|b| ). 

For the subsamples, the early-type galaxies give (3 > 
4.0 with 99.96% confidence (Ax 2 = 11.3 for e = 0.34) 
and (3 > 4.5 with 97% confidence (Ax 2 = 3.8), whereas 
the late-type and power-law galaxy subsamples are each 
consistent with (3 = 4.0 (A% 2 < 1). Core galaxies exceed 
(3 = 4.0 with marginal significance (A% 2 = 1.1) and are 
consistent with (3 = 4.5. Including central kinematics in 
the definition of a further erodes the significance of high 
(3 for the early-type and core galaxy subsamples. More 
robust black hole measurements and more sophisticated 
theoretical models ta king into account of g alaxy types 
and environment (e.g., Zubovas fc King|2012 ) are needed 
before stronger constraints can be obtained. 

The intrinsic scatter in Af. plotted in Figure [5] for 
different intervals of cr, L, and Afbuige serves as an in- 
dependent test for theoretical models of black hole and 
galaxy growth. Our dataset shows decreasing scatter in 
M, with increasing a. However, there are currently in- 
sufficient data to probe the 30— 100 km s^ 1 range, where 
scatter in M, could identify the initial formation mecha- 
nism for massive black holes ( Volonteri et al. 2008 Volon- 
teri fc Natarajan|20 09). Theoretical models of hierarchi- 
cal mergers in ACUM cosmology predict that scatter in 
M, should decline steadily with increasing stellar mass 
(Af*), even when M , and M± are initially uncorrelated 
jMalbon ^ etg]|2Q07| [Peng||2007t prschmann et al ||2010 



Jahnke c? Maccio 2011). The semi-analytic models by 



Malbon et al. (2007), for instance, predict that black 



holes with present-day masses > 10 s Af Q have gained 
most of their mass via black hole-black hole mergers, 
yielding extremely low scatter (eo ~ 0.1) at the upper 
end of the Af, — Afbui ge relation. More recent models by 



Jahnke & Maccio ( 2011 1 use fully decoupled prescriptions 
for star formation and black hole growth, and attain a 
more realistic amount of scatter on average. Yet these 
models still exhibit decreasing scatter as Afbuige increases 
from ~ 10 9 Af Q to ~ lO 115 A/ . In comparison, we ob- 
serve nearly constant scatter from Afbuige ~ 1 11 Mq to 
10 1 2 Af(7), beyond the highes t bulge masses produced in 



the Jahnke & Maccio ( 2011[ ) models. 

Our final comment is that investigations using the 
M m — a correlation should consider the definition of cr, i.e., 
whether it is measured from an inner radius of zero or 
r mi- We find that both definitions yield similar amounts 
of scatter in the M, — a relation (Table [2]) , so neither 
has a clear advantage for predicting Af/ Excluding 
data within r; n f corresponds more closely to cases where 
r; n f is unresolved, such as seeing-limited galaxy surveys, 
high-redshift observations, or numerical simulations with 
limited spatial resolution. From a theoretical perspec- 
tive, the evolutionary origin of an Af. — a relation and 
the immediate effects of gravity may warrant separate 
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consideration. On the other hand, the total gravitational 
potential of a galaxy includes its black hole. Our test of 
how redefining a alters the M. — a relation has only con- 
sidered 12 galaxies for which data within r- m { contribute 
prominently to the spatially integrated velocity disper- 
sion. At present, the full sample of M m and a measure- 
ments comprises a heterogeneous selection of kinematic 
data. Rather than advocating a particular definition, we 
wish to call attention to the nuances of interpreting the 
M m — a relation and encourage future investigators to 
consider their options carefully. 

As this m anuscript was being finalized, |van den Bosch| 
et al. (2012) reported a measurement of M, = 1.7±0.3 x 



nPl^m" NGC 1277. They reported a = 333 km s" 1 
for data between r- ln f and r c g, and Mb u ig e = 1.2 ± 
0.4 x 10 11 M Q ; the latter measurement suggests that 
NGC 1277 lies two orders of magnitude above the mean 
M, -Mbuige relation. Adding NGC 1277 to our 72-galaxy 
sample changes our global power-law fit to the M, — a re- 
lation only slightly: a = 8.33±0.05, (3 = 5.73±0.32, and 
e = 0.39 (from MPFITEXY). Our global fit to M.(M hu i ge ) 
for 36 galaxies including NGC 1277 yields a = 8.51±0.09, 
p = 1.05 ±0.13, and e = 0.44. 

Dynamical measurements of M, require substantial ob- 
servational resources and careful analysis, and are often 
published individually. Nonetheless, recent and ongoing 



efforts are rapidly expanding the available M, measure- 
ments and revising the empirical black hole scaling rela- 
tions. Our online database 4 is aimed to provide all re- 
searchers easy access to frequently updated compilation 
of supermassive black holes with direct dynamical mass 
measurements and their host galaxy properties. Updated 
scaling relations can be used to estimate M, more ac- 
curately in individual galaxies. This can improve our 
knowledge of Eddington rates and spectral energy dis- 
tributions for accreting black holes, as well as time and 
distance scales for tidal disruption events. Moreover, the 
M, — a relation for quiescent black holes has been used 
to normalize the black hole masses obtained from rever- 



beration mapping studies of active gala xies ( Onken et al. 
2004| |Woo et al |2010| |Park et al.|2012| . This important 
calibration could be improved by addressing morphol- 
ogy biases in the reverberation mapping samples and the 
M, — a relations for different galaxy types. 
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Table Al, continued 
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Notes: The first reference column corresponds to the black hole mass measurement, and the second corresponds to the measurement of 
M-t I L used to compute Mb u i ge . Bulge luminosity Ly is in solar units. Quoted errors (+, — ) for M. are 68% confidence intervals. We assume 
0.24 dex uncertainty for all iWb u i go values. The black hole radius of influence rj n f is defined by GAf./cr 2 . Morphologies include designations 
for po wer-law (pi), core (C), an d intermediate (I) surface brightness profiles. Distances for 44 objects have been upd ated since the c ompila- 
tion of |McConnell et al. (2011a); see notes w-z below. A more detailed version of this table is available at http : / /blackhole . berkeley . edu 

References : (1 =| Ghez et al. | 2008[l ; ( 2 =|Gi llessen et al. | 2009]l ; (3 =| Dalla Bonta et al. | 200 9l ; (4 = | GreenhiIl et al. | 2003l ; ( 5 =|Ca ppellari] 
= Vcrolmc et al. 2002) ; (7 = Bender et al. 2005) ; (8 = Krajnovic et al. 2009) ; (9 = Schulzc & Gcbhardt 2011) ; (10 = Bower 
=|Kuoet al.|201lt ; (12 =|Atkmson et al.|20U5ll : (13 =|lNowaket al.|2O08[ ; (14 = Rush et al. 20lT1l ; (15"=jRusli 20121 ; (16 



et al.|2002b ; (6 
et al. 1201)11 ; (11 



= |Gebhardt et al.|2UTj7] ; (17 =lTloughton et al. 2006) : (18 =|Sarzi et al. 2001) ; (19 =|Devereux et al.|2U03|l ; fZO = Emscllcm et al.|1999ll 
; (21 = |TJayies et al.jl2U0b| ; (22 =|Barth et al. l ijOU JL t ~; (23 = |Nowak et al.|201oj ; ( 24 =|van"den Bosch & de Zeeuw|20l0|) ; (25 = |k ondratko| 



et al. 2008) ; (2b = Gultckin et al. 2009b) ; (27 = McComiell et al. 2011a) ; (28 = Walsh et al. 2012) ; (29 = Hcrrnstem et al. 2005) ; (30 : 
1'crrarcsc ct al. 1996) ; (31 = (Jretton & van den Bosch 1999a) ; (32 = Walsh et al. 2010) ; (33 = Gcbhardt et al. 2011) ; (34 = Nowak et al. 



2007) 7735 = | Jardel et al. | 201 l [ ; (36 =|Shen fc Gep hardt 20 lTj]| ; (37 =|Korpqendy et al. | 201l| ; (38 =|de ^'rancesco et al.|2U08|l ; (39 
Cappellari et al.|2UU9^ ; (40 = | McUonnell et al.|2011b| ; (41 = | F£rrarese fc IbrdjiyyCT ; (42 = [van der Marel van den Bosch 1998 ) ; (43 = 
Wold et al.|2U0bj ;(44 =[McConncll ct al.|2012| ; (45 = [Harmg fc Rix|lj004^ ; (46 = |Magorrian et al.|l998| ; (47^ Gcbhardt & Ttiomasj2009]) . 



Notes on individual galaxies: 

a The literature contains a large number of estimates for the velocity dispersion of our G alaxy's bulge, us ng d ifferent kinematic 



tracers at different radii. We use the radially averaged measurement of a = 103 ± 20 kms 1 from |Tremaine et al 



(2002) 



We have re-computed a for 12 galaxies, considering kinematic data between r; n f and r e g. The corresponding values of a are 
listed here. Table [Tj also lists the values of a using data from to r e g. 

c Maser-based black hole mas ses for severa l galax ies are presented in |Greene et al.| ( |2010| l and |Kuo et al.| ( |2011} . We use the 
veloci t y disp ersions presented in |Greene et al. (2010). For consist ency with the rest of our sample, we use the black hole masses from |Kuo| 
|et al.l |2011[l, which agre e with the values in [Greene et al.| (|2010[ l but do not include distance uncertainties in the overall uncertainty tor 
1 1 " provide an updated distance and black hole mass for UGC 3789, which are consistent with the values we adopt 



TT TBraatz et al.| J2010[ > 
>m ]kuo et al.| | |20ll| . 



from|fvuo et al 

d Following G09, our sample includes two distinct measurements for NGC 1399. We weight each of these measurements by 50% 
when performing fits to the black hole scaling relations. 

e The literature con tains two inconsistent estimates of the V-band luminosity of NGC 3607: My = —21.62 in G09, and My = —19.88 in 
|Lauer et aT] ( |2007a| . 

' The literature contains two inco nsistent estimates of the ^/-band bulge luminosity of NGC 4564: My = —19.60 in G09, and 
My = -20.26 in |Lauer etHI] l |2007afr . 

9 Bulge luminosities for NGC 4736 and NGC 4826 were included in the sa mple of [M cConnc ll et al. | ( |2011a| and their fit to the 
M, — L relation. These luminosities corresponded to pseudobulges identified in|Korrncndy ct al. (2011), and we have not included them 
in our present fits. 

h The stellar dynamical measurement of M. in NGC 5128 by |CappelIari et ah] | |2009| is fully consistent with the molecular gas 
measurement M, = 5.3 



+0.6 



(2007) 



w We have updated the di stances to 29 galaxies using surface brightness fluctuation measurements from |Tonry et al.| ( |2001ty , with 
the corrections suggested by Blakeslee et al. (2 010[ ). 

x We have updated the distances to 12 galaxies using surface brightness fluctuation measurements from |Blakeslee et al~] | |2009| , 
which are based on data from ACS on the Hubble Space Telescope. 

y We have adopted a distance of 0.73 Mpc to M31, based on Cepheid variable measurements by |Vilardell et al.] ( |2007[ ) ■ We as- 
sume that M31 and M32 lie at the same distance. 

z For NGC 4342, we have adopted the distance of 23 Mpc by |Bogdan et al.| poT2} . 
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